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ABSTRACT 

Many theoretical and laboratory studies predict H2 to be formed in highly 
excited ro- vibrational states. The consequent relaxation of excited levels via a 
cascade of infrared transitions might be observable in emission from suitable in- 
terstellar regions. In this work, we model H2 formation pumping in standard 
dense clouds, taking into account the H/H2 transition zone, through an accu- 
rate description of chemistry and radiative transfer. The model includes recent 
laboratory data on H2 formation, as well as the effects of the interstellar UV 
field, predicting the populations of gas-phase H2 molecules and their IR emission 
spectra. Calculations suggest that some vibrationally excited states of H2 might 
be detectable towards lines of sight where significant destruction of H2 occurs, 
such as X-ray sources, and provide a possible explanation as to why observational 
attempts resulted in no detections reported to date. 
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Introduction 



By far the most abundant element in the universe is hydrogen. Consequently H2 is the 
most abundant molecule and is the dominant collision partner in dark interstellar clouds. 
Dust grain surfac es act as heterogeneous catalysts in the formation of H2 molecules from 
atomic hydrogen (iGould fc Salpeterlll963l : iHoUenbach &: Salpeterlll979l ) . Beyond this general 
consensus, the actual formation mechanism remains elusive and the internal energy distribu- 
tion of the nascent hydrogen molecule is unknown. During formation, the H2 binding energy 
AEh = 4.476 eV must be partitioned between the ro- vibration al excitation and translat ional 
energy of the nascent molecule and heating of the dust grain (IDuley fc Williams! Il993l ). By 
studying the formation pumping of molecular hydrogen, namely the ro-vibrational distribu- 
tion of nascent H2 molecules, we can constrain interstellar chemistry both in the gas-phase 
and on grain surfaces. 

If the internal and translational energies of nascent molecules are relatively small, then 
significant grain heating must take place, which may lead to the desorption o f volatile 
molecules from the dust grain surface (IDuley fc Williamslll993t [Roberts et al.ll2007| ). The H2 
internal energy distribution could have a significant impact on the chemistry occurring in the 
interstellar medium (ISM) because vibrationally excited H2 will increase the overall energy 
budget of gas-phase processes. There have been many theoretical and l aboratory studies 



that predict to be formed in highly excited ro-vibrational states (e.g.. lSizun et al.ll2010 



Martinazzo fc Tantardinil |2006| and references therein) . It is possible that this formation 
pumping may be observable in the infrared (IR) spectra of H2 molecules. 

The effects induced in the IR spectrum of Ho by ro -vibrational excitation of nascent 
molecules was first considered by lBlack fc Dalgarnd (119761 ) . who employed a formation pump- 
ing model in which equipartition of the H-H binding energy released was arbitrarily assumed. 
In this model, the binding energy is split equally between the internal energy of the molecule, 
its translational energy on desorption from the grain surface and the heat imparted to the 
grain lattice. The molecule is formed at an effective temperature Tf ~ 9000 K. The internal 
energy is spread with a Boltzmann distribution throughout the ro-vibrational levels, with 
the ortho-to-para ratio (OPR) being approximately 3. Subsequently, several classical molec- 
ular dynamics and quantum mechanical calculations have been carried out for H2 formation 
on surfaces whose chemical compositions are analogous to interstellar dust grains, but the 
results have shown a wide dispersion in predicted vibrational distributions. 
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Draine fc Bertoldil (Il996[ ) proposed a ro- vibrational distribution function that boosts 
the populations of vibrational states relative to the population of rotational states. The 
formation temperature is = 5 x 10^ K, the OPR is ~ 2.8 and the mean vibrational and 



rotational levels are 5.3 and 8.7, respectively. iDuley fc Williamsl ( Il986l ) suggested a mech- 
anism which differs from previous proposals in the method of stabilization of the reacting 
complex. The stabilization energy (~ 0.4 eV) is transferred to a surfac e band, whose energy 
is that of the OH stretching vibration. The iDuley fc Williamsl ( 1l986l ) model predicts that 
the H2 molecule on formatio n is ejected into the gas vi brati onally excited (v ~ 6, 7) b ut rota- 
tionally cool (J = . 11. ThelOulev fc Williamsl fll986h and lsiack fc Dakarnol Jl976[) models 



were employed by iLe Bourlot et al.l ( Il995l ) to generate H2 IR spectra. iLe Bourlot et al. 
( 119951 ) presented results for another formation pumping scheme, which postulates that no 
energy is transferred to translation or to dust lattice modes. Thus, H2 is formed in its 
highest vibrational level, v = 14, close to the dissociation threshold, w'ith J = 2 and 3 
weighted by the nuclear spin statistics. H2 spectra were also pres ented bvlTine et ahl (120031) 



based o n ER quantum calcu l ations by iFarebrother et al.l (120001 ) and iMeijer et al. 



(l200lh . 



Finally, iTakahashi fc Ueharal ( I2OOII ) constructed formation pumping models for hydrogen 
molecules newly-formed on icy mantles, carbonaceous and si hcate dust, based on classi- 



cal and quantum theoretical studies of molecular dynam ics (iParneix fc Brechignad Il998 
Takahashi. Masuda. fc Nagaokal Il999l : iMeijer et al.ll200ll ). All these models provide char 



acteristic spectral patterns, which may be used to discriminate between the H2 formation 
pumping mechanisms via astronomical observations. 



Duley &: Williamsl (Il993l ) suggested that the most favourable location to detect forma- 



tion pumping would be from dense, dark, quiescent, sta r-less cores, where ul t raviol et (UV) 



pumping is minim ised. This suggestion is supported by ITakahashi fc Ueharal ( 1200 ll ) and by 



Tine et al.l (120031 ). who state that the relative emissivities of lines due to formation pumping 



in dense clouds can be a factor of 500 greater than in diffuse clouds. However, dark cloud 



Tine et al. 


2003 


and 


Coneiu et al. 


2009) 



sider that non-detections of H2 in dark clouds may be explained by the thermalization of 
nascent H2 molecules on the surface of dust grains, either within the porous structure of the 
grain or via collisions with other adsorbates in the thick icy mantle surrounding dust grains 
in dark clouds. Furthermore, spectra may be not observable due to the lack of an efficient 
mechanism to supply H atoms, such as a very high cosmi c -ray i onization rate ( ~ 10~^^ s~^ 



to dissociate H2 molecules, as found by lLe Bourlot et al.l (119951 ). 



In this paper, we present a new experimentally derived forma tion pumping model 



constructed from the results of the UCL Cosmic Dust Experiment ( iPerry &: Pried 12003 



Creighan. Perry fc Pricdl2006l ) and we compute IR emission spectra of H2 expected in inter- 



-4- 



stellar clouds. The H2 formation pumping model is coupled with an accurate description of 
both radiative transfer and chemistry in stratified dark clouds. Section 2 contains a descrip- 
tion of the formation pumping excitation model employed in the calculations. In Section 3, 
we summarize the method and the procedures followed using different formation pumping 
models to calculate the complete ro-vibrational distribution of H2 formed on dust grains. 
H2 emission spectra are presented in Section 4. We discuss observational implications and 
we present our conclusions in Section 5. 



2. H2 formation pumping 



2.1. The experiment 

The ro-vibrational excitation of molecular hydrogen desorbed from surfaces can be 
measured experimentally. A quantitative partition of the excitation between vibrational 
states has been investigated by the UCL Cosmic Dust Experiment, which probes the dis- 



tribution of the rotational states within each vibrational manifold jPerry fc Pried 12003 



Creighan. Perry fc Pried l2006l ). The experiment studies the formation of molecular hydro- 
gen, primarily HD, on a highly-oriented pyrolitic graphite (HOPG) surface under ultrahigh 
vacuum, following continuous irradiation of the surface by H and D atoms. The nascent HD 
(or H2) molecule will desorb in ro-vibrational states {v, J) of the ground electronic state. 
Hydrogen molecules are state-selectively ionised using laser-induced resonance enhanced 
multi-photon ionisation spectroscopy. The relative populations of the ro-vibrational states 
are then derived from ion yields. 

The experiment probes HD preferentially to H2 (or D2) because there is a significant 
amount of undissociated H2 originating from the H-atom source. As studies of w = 1 and 2 for 
the isotopic species re vealed very similar flux densitie s and rotational distributions for both 
HD and H2 molecules (jCreighan. Perry fc Pricdl2006l ). we extrapolate the HD v = 3 — 7 data 
to obtain the ro-vibrational distribution of nascent H2. However, nascent HD in the f = 
stat e could not be detected above the signal from background gas in the vacuum chamber 
(see ICreighan. Perry &: Pried l2006l ). As a modest estimate of the internal excitation of H2, 
we set the v = populations to be equal to the populations of the v = 1 states. Since this 
is an arbitrary choice, H2 emission spectra have been generated with the v = states both 
significantly more and less populated than expected. The sensitivity of the H2 spectra to 
the error in the population of the v = states is discussed further in Section 4. 

HD formed on HOPG held at 15 K ha s been found in ro-vibrational states (v, J) = (1 — 
2, 0-4), (3-5, 0-6), (6, 0-4) and (7, 0-3) JCreighan. Perrv fc PricJ20od : llslam Latimer fc Price 
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20071 : iLatimer. Islam fc Pried l2008l ). The ro-vibrational distribution of HD peaks ai v = 
4, J = 1. The total vibrational population, found by summing the relative rotational pop- 
ulations, approximately doubles with subsequent state up to f = 4, then falling sharply 
at f = 5. This indicates that the vibrational distribution differs dramatically from a Boltz- 
mann distribution. Examination of Boltzmann plots reveals that the rotational distributions 
within each w— state also slightly deviate from a Boltzmann distri bution. Translational en- 
ergy of the molecule is known to have an upper limit of 0.9 eV ( ICreighan. Perry fc Price 
20061 ). Such a low translatio nal energy of nasc ent hydrogen molecules has been observed by 
other experimental studies (IVidali et al.l 120041 ). From assuming that the f = populations 
are roughly equal to the v = 1 populations, the ro-vibrational distribution corresponds to 
an internal excitation of 1.74 eV, with at least 41% of the HD binding energy flowing into 
the surface. 

In order to retain any deviation from the Boltzmann distribution, we scale the HD 
rotational populations by the appropriate nuclear spin statistical weight Qn obtain the ro- 
tational popujBtioiis_ofJl2^. A formation in space to be approximately 
3 (e.g., iDraine fc Bertoldilll996l and iTakahashil 1200 ll ) . we set = 1 and 3 for even— J and 
odd— J, respectively. By scaling the ro-vibrational populations from HD to H2, we note that 
although the ro-vibrational distribution is preserved, the average energy of the molecule has 
marginally increased. This change in internal energy is due to the energy levels of H2 being 
more spaced out within the potential well, as H2 is lighter than HD. Thus the resulting aver- 
age energy of H2 is 1.95 eV. Three body coupling between the two H atoms and a quantum 
defect within the substrate may be strong enough to randomize the spin orientation of the 
newly formed molecule. However, we do not expect there to be as many defects on our 
HOPG surface as on a real interstellar dust grain (see Sectio n 2.2). Experiments conducted 
probing the ro-vibrational distribution of H2 as well as HD (ICreighan. Perry fc Pried l2006l ) 
show that there is no signiflcant difference in vibrational level populations for H2 and HD 
formed in f = 1 and 2. Also, rotational temperatures of H2 and HD within each vibrational 
state were found to be similar (~ 300 K), taking an OPR equal to 3, for H2. Further- 



more, iTakahashi fc Ueharal (1200 ll ) calculate that the OPR of newly-formed H2 is nearly 3 
on silicate, carbonaceous and icy surfaces. Thus, we feel that statistical weighting is a good 
approximation. 

The relative populations of all the states are then normalised to obtain 6jjj, the fraction 
of H2 formed on grain surfaces that leaves the grain in level {v, J), such that j 6yj = 1. 
The formation pumping population distribution of H2 used throughout this paper is shown 
in Fig. [Hand reported in Table 1. 
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2.2. Astrophysical surfaces 



Astrophysical surfaces are of course likely to be very different from the HOPG surface 
used in the UCL Cosmic Dust Experiment. Important characteristics of an interstellar dust 
grain are chemical composition, fraction of crystallization, roughness, state of charge and 
temperature. Therefore, it is important to discuss how the adopted representation of the H2 
formation pumping population distribution may be related to the properties of interstellar 
grain surfaces. 



The experimental results utilise d in this work (jCreighan. Perry fc Pricell2006l : llslam. Latimer fc Price 



20071 : iLatimer. Islam fc Pried 120081 ) show that HD forms vibrationally e xcited on HOPG, 
which is a non-porous, well-ordered crystalline structure. lAmiaud et al.l (120071 ) have shown 
that D2 forms vibrationally excited [v > 2) on non-porous amorphous ice (NP ASW) held at 
8 — 30 K. Therefore, vibrational excitation of nascent molecules seems to occur regardless o f 
crystallization fraction. Both the UCL Cosmic Dust Experiment and lAmiaud et al.l (120071 ) 
have detected nascent molecular hydrogen in vibrational states v = 1 — 7. The detection 
of vibrationally excited molecules in the two experiments implies that the formation pro- 
cess on HOPG and NP ASW may be similar. Consequently, internal excitation of nascent 
hydrogen molecules may not significantly depend on the chemical composition or degree of 
order of the surface, given a non-porous surface. Unfortunately there is no work to date 
which quantifies the ro-vibrational distribution of nascent molecular hy drogen on astrophys- 
ically relevant surfaces, other than the UCL Cosmic Dust Experiment. iGough et al.l (Il996[ ) 
have shown that H2 forms ro-vibrationally excited on carbon surfaces with temperatures 
of 90 — 300 K, but these temp eratures are too high to simulate the majority of interstellar 
conditions. IVidali et al.l (20061) have found that the recombination of hydrogen is efficient on 
carbon surfaces at temperatures 11 — 18 K, on olivine surfaces at temperatures 6 — 9 K, on low 
density amorphous ice surfaces at temperatures 11 — 15 K and on high density amorphous 
ice surfaces at tempera tures 14 — 18 K. H owever, the temperature-programmed desorption 
(TPD) experiments of IVidali et al.l (120061 ) do not probe the ro-vibrational distribution of 
nascent molecules. 

The roughness of a dust grain may also affect the ro-vibrational distribution of nascent 
molecules, as defects in the crystalline structure may allow atoms to chemisorb, hence bond 
more strongly wit h the graphite su rface. There is a 0.2 eV barrier to chemisorption on 
graphite surfaces (jZecho et al. 2002 ). which arises from carbon atoms having to pucker out 
of the graphite sheet to bond with incident hydrogen atoms. Hence for the UCL Cosmic 
Dust Experiment, where the incident H- and D-atoms are at T ~ 300 K, chemisorption is 
improbable unless the formation mechanism is dominated by reactions at defe cts. TPD of 



an etched graphite surface irradiated by H-atoms at 2000 K was conducted by IZecho et al. 
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( 120021 ) ■ The high temperature of the H-atoms allowed them to overcome the 0.2 eV barrier 
and chemisorb to the surf ace. The e tched surface created terrace edges on the graphite 



lattice, namely defects. IZecho et al.l ( 120021 ) found that atoms bond more strongly at the 
defects. However, there was minimal change in the recombination of chemisorbed atoms with 
the density of terrace edges on the graphite surface, implying that the formation process was 
dominated by conventional sites on the planar surface. Hydrogen atoms adsorbed at terrace 
edges were found to desorb primarily in the form of hydrocarbons, rather than molecular 
hyd rogen. However, in dark clouds, the energy of the incident atoms is much lower than in 
the IZecho et al.l (|2002[ ) experiment. Therefore, it is still feasible that defects dominate the 
reaction of hydrogen atoms at low temperatures. In this cas e, theoretical work based o n 
chemisorption and the Eley-Rideal mechanism, such as that by lTakahashi fc Ueharal ( 120011 ). 
may be more accurate than results from the UCL Cosmic Dust Experiment. 



The charge state of a grain is likely to affect surface chemistry (e.g., ICaruana fc Holt 



20101). If the charge on the grain is largely delocalized, then one would not expect there to 
be much effect on the recombination of two H-atoms. If the charge is localized near the H2 
formation site, then it may be harder for the molecule to escape the grain, hence reducing 
the formation rate. 

A more important factor is the dust grain morphology, in particula r the porosity of 



the surface. Both bare and icy grains may exh ibit a microporous structure (jGreenbergll2002 



Williams &: Herbstll2002l : IWilliams et al.ll2007| ). Alt hough atoms are mobil e on porous grains 
dMatar et al.lboosh and recombination is efficient (IHornekaer et al.ll2003l ). the nascent hy- 
drogen molecules are found to thermalize in the pores, losing kinetic and internal energy. 
Therefore, although the H2 molecule may originally form with ro-vibrational excitation, sub- 
sequent collisions with pore walls may mean that there is no apparent formation pumping 
of molecules that have escaped the grain surface. ICongiu et al.l (120091 ) have shown experi- 
mentally that D2 does not form vibrationally excited on porous surfaces. 

The work presented in this paper incorporates our new formation pumping model from 
the UCL Cosm ic Dust Experiment, as wel l as some theoretical models, such as the work 
put forward by iTakahashi &: Ueharal (I2OOII ). We also include a model in which there is no 
apparent formation pumping. We compare spectra generated by these different formation 
pumping models in order to find observational markers to discriminate between formation 
pumping mechanisms. 
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3. The H2 level distribution model 



We have constructed models to compute the H2 level distribution expected in a ther- 
mally exci ted gas, based on a rad i ative transfer code developed to study H2 formation 
pumpin g JCasu &: Cecchi-Pestellinil 2005 ) and e xcitation in turbulent diffuse interstellar 
clouds ( Cecchi-Pestellini. Casu &: Dalgarnol 20051) . T he ap proach is similar in some respects 
to those developed by Sternberg fc Dalgarnol ( 1989 ) and Draine fc Bertoldi ( 1996). How- 
ever, we do not adopt the cascade efficiency factor formalism (IBlack fc Dalgarndll976l ) and 
we solve the full set of statistical equilibrium equations for the first M = 300 levels of 
H2. We describe the depth-dependent photodissociation r ates by self-shielding func- 
tions ( van Dishoeck fc Black 19881: Sternberg fc Dalgarno 1995), including the prescription 



Draine fc Bertoldi ( Il996[ ). All radiation induced processes have been 



for line overlap given in 
computed taking into account dust extinc tion, which has been assum ed to follow the mean 
galactic interstellar extinction curve (e.g.. iFitzpatrick fc Massal 120071 ) . 



In statistical equilibrium the A/" populations 
to the set of algebraic equations 



{v, J), of the H2 levels are solutions 



n. 



j>i 



rijAji + ^ rij {Cji + Wji) + ICl + TZn^iniSi 



(1) 



In eq.([T]), Wij are the excitation rates from the level i to level j via UV pumping to electron- 
ically excited states, are the Einstein coefficients for spontaneous radiative decay, Cij are 
the temperature-dependent collisional rates, /3j is the rate of photodissociation out of the 
level i, Vi is the rate of additional destruction processes, such as collisional dissociation and 
ionization, and ICf are entry (+) and exit (-) chemical rates. The cosmic-ray destruction 
rate is denoted by C,. The last term in the r.h.s. of eq.([T]) describes the formation of II2 via 
grain catalysis: TZ is the formation rate, tt-h the total volume density of hydrogen, rii the 
volume density of atomic hydrogen, and 5i is the fraction of H2 formed on grain surfaces 
that leaves the grain in level i. The highest lying state is (f , J) = (3, 27) at about 52,000 K 
above the ground state. The level populations are subject to the normalization conditions 
J2vj ''^vJ ~ ''^H2 ; where is the number density of hydrogen molecules. 

The radiative transfer code includes inelastic collisions with H, He and ortho- and 
para-Ho, with fully quantum rn e chani c al calculations of collisional rates given in studies by 
Floweil (ll998h.lFlower fc Rouefj Jl998h. IPlower fc Rouefj (fl998 b[). IPlower. Roueff fc Zeip pen 
(119981 ) and iLe Bourlot. Pineau de Porets &: Plowerl (119991 ) . Por levels where quantum calcu- 
lations are not available, the extrapolation scheme for the H-H2 rate collisions provided 
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by Le Bourlot and co-workers in their code for photon— dominated regions (PDRs) has 
been adopted (Le Bourlot private communication). For all other coUisio nal partner s when 
quantum calculations are lacking, the collision scheme put forward by iTine et al.l (Il997l ) 
has been employed. For models with temperatures larger than 30 K, we have i n corpo - 
rated three rate coefficients for ortho-para conversions provided by lSun fc Dalgarnd ( 1l994j ). 
Energy levels, radiative decay rates and diss oc iation probab il ities f o r electronic transitions 
hav e been published by kbgrall et al.l Jl992[ l. kbgrall et all Jl993h . kbgrall et all Jl993bl ) 
and kberall. Roueff fc Driral J200oh . Extra data, covering levels up to J = 25, were kindly 

provided by Abgrall (private communication). Quadrupole radia tive decays and energies of 

ro-vib rational levels of the ground electronic state were taken from lWolniewicz. Simbotin fc Dalgarno 
fll998h . 



The H2 formation rate is more fully described by moment equations (iLe Petit et al. 



2OO9I ). rather than by rate equations using the expression 7lnu_ni6i. This is because small 
grain sizes and low atom fluxes are subject to large flu ctuations, and t hus c alculating the 
H2 formation rate requires stochastic methods. However, iLe Petit et al.l (120091 ) find that the 
moment equation results agree with the rate equation results in a wide range of conditions, 
except for dust grains at temperatures larger than 18 K, in which case the rate equations 
overestimate the H2 formation rate. For dark clouds, as investigated in this paper, where 
dust temperature are ~ 10 K throughout the cloud, the H2 formation rate can be ade- 
quately described using the standard rate equation term assumed in eq.([T]). Of course, if the 
radiation environment differs substantially from the standard interstellar conditions, dust 
temperatures may be noticebly larger than 20 K (see Section 4). 

In the present study, for all formation pumping models, except the one presented in 
Section 2, the relative (normalized) populations of the ro-vibrational states are given by 

= Cgr,fi{v,J)exp[-f2{AE,j,Tf)] (2) 

where C is a normalization constant and AE^j is the energy in K of level {v, J) referred to 
the ground state, /i and /2 are shape functions depending on the specific formation model. 
We consider the following H2 formation pumping models: 



(i) the formation pumping model described in Section 2, where 6yj is taken directly from 
extrapolation of HD experiments (see Table 1); 



the acquired internal energy, Ei is statistically distributed among the energy levels, 
/i = 2J + 1 and /2 = A E„j/Tf; by setting the fo rmation temperature to Ty = 9,000 K 
the model proposed by iBlack fc Dalgarnd ( 1l976l ) is recovered 

X (2 J + 1) AE^j exp {-AE,j/Tf) 1.5 eV 



Qn X (2 J + 1) exp {-AE,j/Tf) 



Kb 



(3) 
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Kb being the Boltzmann constant; 



{Hi) to enhance the populations of high v states with respect to high J states, we set 
/i = f + 1; this clas s of formation models wi th T/ = 50, 000 K provides the pumping 



profile suggested by iDraine fc Bertoldil ( 119961 ) ; 



[iv) the three ro- vibrational population distributions for H2 newly-for med on carbonaceous 
dust { iv — c), sihcate dust {iv — s) and icy mantles {iv — i) given bv|Takahashi fc Uehara 
(2001), model A; the functions /i and /2 are taken from Takahashi fc Uehara fcoOl ): 



the same as in {iv) but for iTakahashi fc Uehara] ( 120011 ) model B; 



[vi) a minimal H2 formation pumping model which limits the ro- vibrational distribution to 
the lowest possible levels (f, J) = (0,0) and (0,1) in the ratio 1:3; /i = /2 = except 
for (f , J) = (0, 0) and (0, 1) when fi = 1; this represents the "no formation pumping" 
case. If nascent H2 molecules in dark clouds ther malize with the dust grain pores or 
the icy mantle, as put forward ICongiu et al.l (|2009[ ). then no formation pumping would 
be detected. 



Ho IR emission from dark clouds 



The cloud is assumed to be a plane-parallel slab of constant density, nn- The cloud 
is two-sided illuminated by an isotropically incident UV radiation field. The fi eld inten- 
sity is assumed to be the UV field scaling parameter, y, times the Drainel (I19781) es timate 
of the mean interstellar radiation field as reported in Sternberg fc Dalgarnol ( 1989 ). The 
attenuation of the field due to dust extinction is computed using the analytical solution 
to the transport equation in plane-parallel geometry given in iFlannery. Roberge fc Rybicki 
( ll980l ). The H/H2 ratio calculations as a function of depth have been su pplemented wit h 
the time- and depth-dependent gas-gr ain UCL_CHEM chemical model (IViti et al.l 120041 ). 
In addition, we use the UCL_PDR code (IBell et al.ll2006l ) in order to get the cloud thermal 
profile. We find that gas temperatures are higher than ~ 10 K only at the very edge of dark 
clouds. However, in the case of translucent clouds, or non-standard illumination (x > 1), 
gas temperatures may be significantly larger than 10 K throughout the cloud. 



We perform calculations for translucent and dark clouds of number densities nn 



10^ — 10^ cm ^, cosmic-ray ionization rate ( 



3 X 10 s ^ and column densities A^h 



4.8 X 10 — 1.6 X 10 cm , c orresponding to visual integrated magnitudes Ay = 3 — 
Snow fc McCall 20061). Dark clouds were chosen following the suggestion by 



100 



Duley fc William^ ( 19931 ) and Tine et al. (2003) that quiescent, dense clouds with no UV 
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pumping are the regions of the ISM best suited to detect formation pumping. We extend 
the calculations to clouds with moderate extinctions in order to get larger formation pump- 
ing rates, although contamination from radiative pumping may be substantial (as is actually 
the case). For such clouds, the UCL_PDR code in the density range 10^ — 10^ cm~^ provides 
thermal profiles in which the "10 K edge" is reached at Ay = 1 — 2.5 mag, measured from 
the external boundary of the cloud. For x ^ 100, the UCL_PDR code also finds that dust 
temperatures are larg er than 20 K. In tha t case the use of rate equations for H2 formation 
may be not accurate (ILe Petit et al.ll2009l ). 

In Fig. [2], we compare the results for the six formation pumping profiles, described 
in Section [3l for the case of dark clouds of hydrogen density nn = 10® cm~^, 7^ = 3 x 
10^^^ cm'^ s~^, C = 3 X 10^^^ s~^, standard UV field (x = 1) and integrated visual extinction 
Ay = 100 mag, hereafter called the reference cloud model (RCM). The emissivities roughly 
scale with the simulated instrumental resolution. We assume a FWHM of 5 x 10^^ fim, 
approximately corresponding to the instrument resolution of echelle spectrometers used by 
current ground based telescopes such as the Phoenix instrument on the Gem ini telescope, 



which has a resolution of 50000 — 80000 for near-IR wavelengths 1 — 5 ( iHinkle et al. 



2OOOI ). All of these pumping models provide distinct spectral features. Therefore, in principle. 



astronomical observations could be used to identify the actual formation pumping mechanism 
taking place in the ISM. However, a common pattern is apparent in all the generated spectra 
in Fig. O As the adopted pumping mechanisms result in quite different H2 level excitation 
rates, the underlying common level distribution must be generated by some global ambient 
mechanism, such as thermal pumping. However, coUisional excitation cannot be significant as 
the RCM kinetic temperature is approximately 10 K in most parts of the cloud. These spectra 
are the result of radiative transfer throughout the whole cloud, enfolding both translucent 
and dense cloud regimes. Therefore, the features present in the IR spectra might originate 
via UV pumping close to the edge of the cloud, in a region in which all the hydrogen is 
not yet in molecular form. In the bottom right panel of Fig. [21 we show the IR spectrum 
generated by model {vi), in which formation pumping is not active, but nevertheless a rich 
spectrum is present: radiative pumping in the translucent regime tends to dominate the 
emission even at large extinctions. This interpretation is supported by the spectrum shown 
in Fig. [31 derived for a cloud model with the same physical parameters as the RCM, but 
with visual extinction Ay = 5 mag. The increase of the emission in the lower Ay cloud 
model with respect to the RCM is mainly due to residual UV radiation coming from the 
opposite edge of the cloud. In Table 2, we report wavelength-integrated emission intensities 
produced in different cloud models, including the spectra displayed in Figs. [2]and[31 Since the 
residual field scales with the size of the cloud (in plane-parallel clouds of constant density), 
an increase in Ay produces a decrease in the emitted integrated intensity, when the main 
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pumping mechanism is radiative. We also note that in the RCM, depending on the formation 
model, H2 formation pumping produces an IR excess of about 20 — 40 % with respect to the 
UV-pumped IR background. 

To highlight spectral features arising from specific H2 formation pumping profiles we 
show in Fig. H] the residual spectra computed for the RCM after subtraction of the "back- 
ground" UV pumping contribution, i.e. the spectrum arising in model (vi) (right bottom 
panel of Fig. [2]). Since radiative transfer couples different parts of the cloud, residual spectra 
may provide an indication of the formation pumping effect, but in general they cannot be 
considered as an "exact" measure of the formation pumping contribution to line excitation. 
The most intense residual emission lines produced by the formation pumping are reported 
for each model in Table 3. Model [i) shares some features with model [in) only. All the 
other formation pumping present several common features. Our proposed formation pump- 
ing model produces spectra where the highest vibrational le vel is = 4, with no line s from 
high rotational states in the 4 — 5 /im region. In model {ii) (IBlack &: Dalgarndll976l ). high 
rotational levels are pumped, as expected, with the highest being J = 16 (at an intensity 
level of 1 X 10""^ ergs s~^ cm~^ sr~^ /im~^). A lthough p resenting the highest int egrated 



intensities, models (ii) (IBlack fc Dalgarndll976l ) and (Hi) (jPraine fc Bertoldil 119961 ) do not 
show very prominent spectral features, because of the wide disp ersion of the internal energy 
over a large number of vibrational states. Models (iv) and (v) (ITakahashi &: Ueharall200ll ) 
show high rotational states, as well as a combination of moderately high v and J, such as 
the (f , J) = (3, 9) state. These models provide quite different spectral patterns, with the 
former showing brighter lines, while the latter exhibits a richer spectrum. The most intense 
transition for model (i) is (4 — 2) 0(3), whereas the (1 — 0) S(7) line is strongest for all other 
models. In Table 3, we present emission lines computed for model (vi). Since, in this model, 
formation pumping is suppressed, spectral features arise from radiative and coUisional pump- 
ing close to the edge of the cloud. Transitions involving low rotational states (J < 5) show 
a systematic mixing of internal and environmental pumping mechanisms. As a consequence, 
our formation pumping model, which predicts little rotational excitation, appears to produce 
an emission spectrum contaminated by external factors. This contamination does not occur 
for cases where the emissio n is dominated by transit ions from high rotational states, such as 
the (1 - 0) S(7) line in the iBlack fc Dakarnol Jl976h model. 



In Figs. [5]and[6l we show the volume emissivities of the (4 — 2) 0(3) and (1 — 0) S(7) 
lines as functions of the optical thickness within the cloud. The (4 — 2) 0(3) transition has 
been computed for models (i) and (vi) in the translucent and dense cloud regimes, while 
the (1 — 0) S(7) transition has been computed for models (ii) to (v). In the figures, we also 
present the most important relative contributions to the population of the upper state of 
the transitions due to the included excitation mechanisms, see eq.([T]); these mechanisms are 
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UV pumping by fluorescent cascade from excited electronic states, IR cascade within the 
ground state, thermal collisions and excitation due to H2 formation on dust grains. The 
excitation of the (4 — 2) 0(3) line is driven by radiation close to the edge of the cloud, 
but is dominated by formation pumping in the cloud interior. Since collisional rates scale 
with density, thermal collisions contribute to the line excitation in response to the very high 
density adopted for the RCM. In the case of the low-density translucent cloud model, thermal 
excitation is negligible. In model (vi), line emission is essentially the same as for model (i) 
until the external UV radiation density declines sharply for Ay > 1 mag. No line excitation 
is produced because formation pumping is suppressed in model (vi). For the transition 
(1 — 0) S(7) in models (m), (iii), (iv) and (v), the situation is different. The line excitation 
at large visual depth is almost completely dominated by the IR cascade (~ 90 %) rather than 
by direct excitation from the formation process (~ 10 %). However, since UV pumping is not 
effective for large Ay, the high energy levels must be populated by formation pumping. The 
IR cascade then populates {v, J) = (1,9), and lower energy levels. This is directly dependent 
on the distribution of the internal energy over a large number of ro-vibrational levels, as 
occurs in models {ii) to (u). In contrast, only a small number of ro-vibrational levels are 
populated in model (i). Therefore, the excitation of the (1 — 0) S(7) hue is essentially driven 
by formation pumping. This conclusion is supported by the lack of significant emission at 
the corresponding wavelength in model (vi) (see Table 3). 

We now discuss the dependency of the IR spectra on cloud parameters. In Table 4, we 
present fine intensities for 4 — 2 0(3), 1 — S(7) and 1 — Q(l) transitions using the set 
of cloud models defined in Table 2. The physical conditions of such models cover a wide 
portion of the parameter space. In Table 5, we show the excitation and de-excitation rates 
for the upper levels of the 4 — 2 0(3) and 1-0 S(7) transitions, for selected positions within 
the cloud. We find that 

• the integrated spectrum intensity increases almost linearly with x, as expected if the 
excitation is dominated by the emission in the translucent regime (see Table 2); the IR 
excess provided by formation pumping decreases with increasing UV radiation density, 
e.g., in the case of model (i) the excess is 20% when x = 1, while it is less than 4% 
when X = 1000; 

• integrated intensities decrease with increasing Ay (and hydrogen column density): 
this is a direct consequence of the fact that, in a double-sided illuminated cloud, 
the boundaries may be reached by residual radiation coming from the opposite edge; 
integrated intensities do not decrease when the gas density is = 10^ cm"^: in such 
a case, since only a partial conversion of hydrogen from atomic to molecular form 
takes place, lines pumped by formation get significantly brighter, even if the overall IR 
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emission does not increase; 

• line spectra vary marginally with increasing number and column densities as soon as 
hydrogen is almost totally converted to molecular form; in clouds where partial H/H2 
conversion occurs, line intensities are larger than in dense clouds where emission paths 
are considerably longer; 

• the strength of a line increases steadily, although not linearly, with kinetic temperature, 
because as the temperature rises, more colhsional pumping occurs, although the shape 
of the spectrum does not greatly alter; in general, when kinetic temperature increases 
over 100 K in a substantial fraction of the cloud, there is a sharp rise in the strength 
of the 1 — Q(l) line, due to thermal excitation of H2 to the (f , J) = (1, 1) level; 

• for both levels {v, J) = (1, 9) and (4, 1) in any position within the RCM, the major 
exit channel is radiative de-excitation; the major entry channel is direct formation 
pumping for level (f , J) = (4, 1) and IR cascade from upper levels for {v, J) = (1, 9); 
since both UV fluorescence and thermal collisions are neglible in most parts of the 
cloud, the pumping of high ro-vibrational states is due H2 formation on dust grains; 

• as there is uncertainty in the population of the v = levels (see Section 2.1), H2 spectra 
have been generated for a distribution where the v = levels are twice as populated as 
the V = 1 levels and for a distribution where the v = levels are half as populated as 
the V = 1 levels; these changes make only a minimal difference to line strengths, with 
spectra similar to within 4%. 

Finally, in Fig. [TJ we show the residual emission spectrum of a standard translucent 
cloud (riH = 10^ cm~^. Ay = 5 mag, all other parameters are as for the RCM) for the case 
of model (i). It appears that in a cloud of moderate density and extinction, the intensity of 
the emitted spectrum is 3 — 4 times more intense than in the case of the RCM, in which the 
hydrogen column density is 20 times larger. The increase in line strength is more evident 
for lines produced during H2 formation. This reflects the much larger abundance of atomic 
hydrogen (by about a factor of 1000) in the lower density cloud. The comparison with the 
RCM case also shows that the very long emission path inside a dense cloud produces only 
a marginal increase in the line intensities, since the abundances of the upper levels in the 
emitting transitions fall abruptly as soon as H2 formation saturates. 
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5. Discussion and Conclusions 

In this work we investigate tlie effects of formation pumping in IR H2 emission spec- 
tra with a new formation pumping model (Table 1). We construct radiative tranfer and 
chemical models for H2 newly-formed on dust grains. By using realistic space-dependent 
cloud models, we find that UV radiative pumping dominates the emission even in clouds 
with very high visual extinction. After eliminating radiative pumping, we obtain residual 
IR spectra due to formation pumping. When cloud kinetic temperatures rise over 100 K, 
thermally excited ro-vibrational levels of the v = 1 manifold may contribute to the overall 
IR spectrum. 

Surprisingly, spectra show a very modest increase with both volume and column den- 
sities. IR emission is expected to scale with the square of volume density and linearly with 
column density, via the formation rate and the path along the line of sight. However, deep 
within a cloud {Ay ^ 2.5 mag) an almost total conversion of hydrogen from atomic to 
molecular form occurs. Thus, the emissivity produced in the inner zone of a cloud is "lost" 
within the contribution originating in the transition zone, in which H2 abundances are still 
comparable to those of atomic hydrogen. In other words, the bulk of the emission is coming 
from the "translucent" regime close to the edge of a cloud. 

This can be easily understood considering a simple two-level model representing an 
emission line u — )■ / arising during the H2 formation process within a homogeneous cloud. In 
such a representation, the line emissivity is given by 

he 

Auinu4>\ with Auifiu TZrinriiSu, (4) 



47rA 



(see Table 5), while the brightness is obtained by integration of eq.(jl]) along the line of sight 
AL, 

he 

six ~ -^;^Tln^xi5u(l)x^L. (5) 

In eqs.(jl]) and ([5]), (px is the assumed line profile, riu is the population of the upper level of 
the transition and Xi is the fractional abundance of atomic hydrogen. The intensity ratio 
between emissions originating from dense and translucent regions in the cloud then results 

We get 5/f ~ 51^ when N§ ~ 1000 x A^^;^. Since in the RCM, N^/N^ ~ 20, we finally obtain 
61^ < 2% 61^. Thus, emission from the translucent outer regions of a cloud dominates the 
spectral line intensity, with only a small contribution from the dense central regions of the 
dark cloud. Consequently, models of formation pumping that are rotationally cool, such as 
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the one put forward in this work, are heavily contaminated by pumping processes arising 
from nebular physical conditions. 

It is possible that only lines from high rotational levels can be identified observationally 
as being due to formation pumping. Unfortunately, most of the models proposed in litera- 
ture have common features (see Table 3) and, in general, it appears difficult to discriminate 
between them. In contrast, our new formation pumping model produces a spectrum differing 
significantly from the line patterns of all the other proposed models, and may provide an 
unambiguous signature for detection, but only from regions not contaminated by UV radi- 
ation. In addition, the observation of H2 molecules newly-formed on dust grains appears 
to be currently very difficult. The peak intensity of the spectral lines are of the order of 
0.001 erg s~^ cm~^ sr~^ fim~^ (see Fig. 2), or ~ 1 x 10~^® W m~^ arcsec"^ /im~^. Using 
a resolution R = A/ A A = 37000 (taken from UKIRT webpage^ we obtain an integrated 
brightness at 2 fim of / ~ 1.2 x 10^^^ W m~^ arcsec"^. Hence, using the UKIRT 1-pixel 
slit, which is 0.609 arcsec wide, we obtain a flux of F ~ 4.5 x 10~^^ W m~^ (for an ex- 
tended source). This flux is much lower in comparison to the 3a 30-minute sensitivity per 
pixel obtained with the echelle grating, which has values 8 x 10~^° W m~^ at 1.6 fim and 
6 X 10~^° W m~^ at 2.2 /im. As a consequence, dark clouds are not a good place to look for 
signatures of H2 formation, since an increase in the emission path does not correspond to a 
signiflcant increase in the intensity of the spectrum. 

Therefore, II2 formation pumping may be undetectable in those regions in which molec- 
ular hydrogen is not destroyed at a fairly fast rate. Thus, within the current regime of 
instrument sensit ivity, the non-detections of H2 IR emission in dense clouds reported by 



Tine et al.l (120031 ) and ICongiu et al.l (120091 ) do not need to be explained by t he thermaliza 



tion o f nascent H2 molecules on the surface of dust grains, as suggested by ICongiu et al. 



(120091 ). In dense PDRs, UV radiation maintains a substantial level of atomic hydrogen in 
the gas. However, the radiation also causes excitation of H2 to the Lyman and Werner bands, 
inducing IR fluorescence to further complicate the emission process. Much more suitable re- 
gions for observing formation pumping appear to be X-ray dominated regions (XDRs). High 
energy X-rays penetrate much deeper into gas clouds than UV photons. X-rays are pref- 
erentially absorbed by heavy elements to produce multiply-charged ions and photo-ionize 
the gas deeply within a cloud. High energy primary photo-electrons deposit their energy 
into the gas, inducing a secondary electron cascade, which ionizes, excites and dissociates 
atomic and molecular species, and also heats the gas through Coulomb collisions. Moreover, 
although the electrons efficiently destroy molecular hydrogen, electron discrete interactions 



1 



http: / /www .jach.hawaii.edu/UKIRT / instruments / cgs4 / optical/resolution. html 
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only provide excitations in the H2 vibrational ladder up to f = 2 (jPalgarno. Yan fc Liulll999l ). 
The radiative decays of the electronically excited H and H2 produces FUV photons, H2 Ly- 
man-Werner photons and H Lya photons. We can estimate the impact of photo-electron 
induc ed UV radiation i n the following way: assuming photo-electrons of mean energy 30 eV 



as m 



Tine et al.l 119971 ) . the number of excitations to the states -B^S+ and C^II", including 



the contrib ution of cascading from hig her singlet states, are approximately 0.6 and 0.4, re- 
spectively ( jPalgarno. Yan fc Liul 119991 ). The excitation rate to the excited electronic states 
is thus ~ (x, where (x is the total X-ray ionizati on rate. The radiative e xcitation of Ly- 



man and Werner bands is roughly 5 x 10 s ^ (jPraine &: Bertoldilll996[ ). We obtain the 



equivalent UV field scaling factor by means of the relation x = Cx/5 x 10^^°. Thus, the 
most prominent feature in the IR spectrum generated using our formation pumping model, 
namely the transition 4 — 2 0(3), would be minimally affected by X-rays, as long as the 
ionization rate is ^ 5 x 10"^" s~^. 

In Fig. IHl we report the intensity of the line 4 — 2 0(3) as a function of the H2 ionization 
rate in the RCM. It is evident that the transition 4 — 2 0(3) might be ob s erved if the H2 
ionization rate is C ^ 1 x 10~^^ s~^, in agreement with iLe Bourlot et al.l ( 1l995l ). In this 
case, since the line intensity scales approximately with the ionization rate, the line would 
be largely dominated by formation pumping rather than radiative pumping. The crucial 
point is that XDRs may mantain large ionization rates up to hydrogen column densities o f 



A^H ~ 10 — 10 cm for X-ray energies larger than 1 keV (ICecchi-Pestellini et al.ll2009l ). 



In conclusion, we model H2 formation pumping in standard dense clouds, taking into 
account the H/H2 transition zone. The model, which includes recent laboratory data on 
H2 formation, as well as the effects of the interstellar UV field, predicts the populations 
of gas-phase H2 molecules and their IR emission spectra. Calculations suggest that some 
vibrationally excited states of H2 might be detectable towards lines of sight where signifi- 
cant destruction of H2 occurs, such as X-ray sources. These results also provide a possible 
explanation of the lack of detection to date of II2 formation pumping in dark clouds. 
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Table 1: Normalized population distribution 
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("^rotational populations oi v — Q are taken equal to the populations of ?; = 1 (see text). 




Fig. 1. — The new formation pumping model. The relative ro-vibrational populations for H2 
are extrapolated from experimental studies of HD (see text) by assuming an ortho— to— para 
ratio of 3 and choosing 5qj = 6i^j. 
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Table 2: Wavelength-integrated 
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Table 3: Strongest emission lines 
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Z.4Uoz4 


(1 - 0) g(i) 


A V 

4. / 


1.9 






ii.2 
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(z — Ij U(zj 


1 

1.3 








14.9 


z.oUzz/ 


(2 - 1) Q(ll) 






i.o 






2 97273 


fo - l"! 0('3'1 


6 2 


2.0 






12.6 


3.54591 


(0 - 0) S(16) 




1.2 








3.72256 


(0 - 0) S(14) 




2.2 








3.83849 


(1 - 1) S(15) 




1.2 








3.84413 


(0-0) S(13) 




1.2 








4.07916 


(2-2) S(15) 




1.3 








4.17881 


(0-0) S(ll) 




2.4 


2.1 


1.3 




4.22152 


(1 - 1) S(12) 




1.2 


1.2 






4.41439 


(1-1) S(ll) 




2.2 1.1 


3.0 


1.6 




4.95163 


(1 - 1) S(9) 




3.3 1.9 


5.6 


2.6 





("^peak line intensities larger than 1 x 10 ergs s ^ cm ^ sr ^ fjum 
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Table 4: Peak line emissions for the cloud models reported in Table 2 
cloud model^"^ intensity 



(10 ^ ergs s ■'^ cm ^ sr -"^ //m ^) 







(0 






(ii) 






(iii) 








L2 


L3 


LI 


L2 


L3 


LI 


L2 


L3 


1 


19.9 


0.05 


16.0 


12.7 


6.5 


11.7 


14.4 


4.6 


13.1 


2 


274 


17.1 


354 


227 


54.9 


327 


239 


42.9 


336 


3 


3407 


204 


4082 


2880 


537 


3777 


3024 


420 


3886 


4 


19977 


32133 


124859 


15834 


33556 


123270 


17034 


32984 


123902 


5 


54.6 


0.19 


38.8 


23.9 


30.2 


22.2 


30.4 


21.1 


27.3 


6 


71.6 


0.17 


46.0 


17.0 


53.4 


16.6 


28.4 


37.2 


25.6 


7 


24.4 


0.11 


19.8 


19.3 


5.1 


17.1 


20.4 


3.6 


17.9 


8 


16.4 


0.07 


13.2 


12.8 


3.6 


11.3 


13.6 


2.5 


11.9 


9 


22.0 


0.08 


17.7 


17.2 


4.7 


15.1 


18.1 


3.3 


15.9 


10 


14.8 


0.05 


11.9 


11.4 


3.4 


10.0 


12.1 


2.4 


10.6 


11 


24.0 


0.07 


20.2 


18.9 


4.6 


17.2 


20.2 


3.2 


18.2 


12 


16.2 


0.04 


13.7 


12.6 


3.3 


11.6 


13.5 


2.3 


12.2 






(iv — c) 






(v — c) 






(vi) 






LI 


L2 


L3 


LI 


L2 


L3 


LI 


L2 


L3 


1 


12.5 


10.7 


11.2 


12.7 


6.1 


12.1 


12.5 


0.04 


11.3 


2 


225 


79.2 


324 


226 


52.5 


329 


225 


16.4 


322 


3 


2855 


752 


3742 


2877 


520 


3816 


2864 


193 


3725 


4 


15779 


34549 


123085 


15975 


33603 


123764 


15968 


31752 


122916 


5 


23.1 


49.9 


20.8 


23.7 


28.3 


23.3 


22.9 


0.14 


20.5 


6 


15.6 


88.2 


14.0 


16.7 


49.9 


18.6 


15.5 


0.09 


13.8 


7 


19.2 


8.4 


16.9 


19.3 


4.8 


17.3 


19.1 


0.10 


16.8 


8 


12.7 


5.9 


11.2 


12.8 


3.4 


11.5 


12.6 


0.06 


11.1 


9 


17.0 


7.8 


14.8 


17.1 


4.4 


15.2 


17.0 


0.07 


14.8 


10 


11.3 


5.5 


9.8 


11.4 


3.2 


10.1 


11.3 


0.05 


9.8 


11 


18.8 


7.6 


16.9 


18.9 


4.3 


17.4 


18.8 


0.06 


16.9 


12 


12.5 


5.5 


11.3 


12.6 


3.1 


11.7 


12.5 


0.04 


11.3 



Numbers refer to Table 2 
Wl: (4 - 2) 0(3), L2: (1 - 0) S(7), L3: (1 - 0) Q(l) 
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Table 5: Entry and exit rates of the upper levels of the transitions (4—2) 0(3) and (1 — 0) S(7) 

computed for the RCM 

rates (s~^) 



Av 












{v,J) 














(mag) 




(4,1) 








(1,9) 












( 


i) 






(ii) 


(Hi) 


{iv 


-c) 


{v 


-c) 








entry 


exit 


entry 


entry 


entry 


entry 


exit 


0.1 


1.3( 


-16) 


3.0( 


-17) 


2.6(-19) 


1.9(-18) 


1.2( 


-17) 


2.0 


rl7) 


3.0( 


-17) 




6.1( 


-16) 


3.1( 


-06) 


1.7(-16) 


1.3(-16) 


3.0( 


-16) 


1.4 


;-16) 


7.1( 


-07) 




3.8( 


-16) 


1.5( 


-09) 


1.5(-19) 


1.0(-19) 


1.9( 


-19) 


1.2 


'-19) 


1.2( 


-09) 




1.2( 


-16) 


2.0( 


-15) 


9.9(-18) 


7.1(-18) 


1.7( 


-17) 


9.5 


;-18) 


3.7( 


-08) 


0.5 


6.4( 


-19) 


3.0( 


-17) 


1.1(-19) 


1.1(-19) 


6.4( 


-19) 


1.1 


;-i8) 


3.0( 


-17) 




2.6( 


-17) 


3.1( 


-06) 


9.5(-18) 


7.3(-18) 


1.7( 


-17) 


7.8 


;-i8) 


7.1( 


-07) 




1.6( 


-17) 


2.5( 


-10) 


9.5(-22) 


6.9(-22) 


1.2( 


-21) 


7.5( 


'-22) 


1.5( 


-10) 




5.3( 


-18) 


1.1( 


-16) 


3.8(-20) 


2.6(-20) 


1.2( 


-19) 


3.5 


'-20) 


2.9( 


-09) 


1.0 


3.0( 


-18) 


3.0( 


-17) 


2.2(-17) 


6.7(-20) 


2.7( 


-19) 


4.6 


;-i9) 


3.0( 


-17) 




1.9( 


-18) 


3.1( 


-06) 


3.9(-18) 


3.0(-18) 


6.9( 


-18) 


3.2 


;-i8) 


7.1( 


-07) 




7.4( 


-19) 


3.5( 


-11) 


3.3(-23) 


2.8)-23) 


4.3( 


-23) 


2.5 


'-23) 


1.5( 


-11) 




7.1( 


-19) 


4.3( 


-17) 


8.2(-21) 


5.3(-21) 


1.3( 


-20) 


7.3 


'-21) 


1.7( 


-09) 


2.5 


3.0( 


-18) 


3.0( 


-17) 


5.6(-21) 


3.4(-20) 


2.5( 


-19) 


4.3 


;-i9) 


3.0( 


-17) 




2.5( 


-22) 


3.1( 


-06) 


3.7(-18) 


2.8(-18) 


6.5( 


-18) 


3.0 


;-i8) 


7.1( 


-07) 




2.8( 


-22) 


1.6( 


-13) 


5.5(-26) 


9.6(-26) 


5.6( 


-26) 


3.1 


'-26) 


2.4( 


-14) 




4.7( 


-19) 


4.7( 


-17) 


3.3(-21) 


5.2(-21) 


1.2( 


-20) 


6.5 


r21) 


1.6( 


-09) 


5.0 


3.0( 


-18) 


3.0( 


-17) 


5.6(-21) 


3.4(-20) 


2.5( 


-19) 


4.3 


;-19) 


3.0( 


-17) 




3.5( 


-19) 


3.1( 


-06) 


3.7(-18) 


2.8(-18) 


6.5( 


-18) 


3.0 


;-i8) 


7.1( 


-07) 




2.3( 


-27) 


3.9( 


-17) 


1.5(-29) 


5.2(-29) 


4.3( 


-30) 


2.2 


;-3o) 


1.0( 


-18) 




4.7( 


-19) 


4.7( 


-17) 


3.3(-21) 


5.2(-21) 


1.2( 


-20) 


6.5( 


;-2i) 


1.6( 


-09) 


10. 


3.0( 


-18) 


3.0( 


-17) 


5.6(-21) 


3.4(-20) 


2.5( 


-19) 


4.3 


;-i9) 


3.0( 


-17) 




3.5( 


-19) 


3.1( 


-06) 


3.7(-18) 


5.2(-29) 


6.5( 


-18) 


3.0 


;-i8) 


7.1( 


-07) 








4.6 


;-24) 














4.6( 


-27) 




4.5( 


-19) 


3.9( 


-17) 


3.3(-21) 


5.2(-21) 


1.2( 


-20) 


6.5 
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1.6( 


-09) 


50. 


3.0( 
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3.0 


;-i7) 


5.6(-21) 
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2.5 


;-19) 


4.3 


;-i9) 


3.0( 


-17) 




3.5( 


-19) 


3.1( 


;-06) 


3.7(-18) 


5.2(-29) 


6.5( 


;-i8) 


3.0( 


;-i8) 


7.1( 


-07) 



T.j>iAjinj/n2 



4.5(-19) 3.9(-17) 3.3(-21) 5.2(-21) 1.2(-20) 6.5(-21) 1.6(-09) 
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Fig. 2. — IR emission in ergs s ^ cm ^ sr ^ /xm ^ computed for the RCM and six dif- 
ferent formation pumpi ng m odels: (i) the model pro posed i n this work (Table 1); (ii) 



lerent lormation pumpi ng m odels: (t) the model pro posed i n tnis worR ( lable 1 ; (ttj 
Black fc Dal garnolJl976h: fmlbraine fc BertoldilJl996h : fw-cl lXakahashi fc Ueh^j200lf ) 
- carbon model A; (f — c) Takahashi &: Ueharal (2001) - carbon model B; (vi) no formation 
pumping. In the case of the iTakahashi &: Ueharal (120011 ) models, we show only the spectra 
for carbon materials, since all the spectra appear to be very similar regardless of the chemical 
composition of grain surfaces. 
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Fig. 3. — IR emission in ergs s ^ cm ^ sr ^ /xm ^ with formation pumping model (i), 
computed for the RCM but with visual extinction reduced to Ay — 5 mag. 
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Fig. 4. — Residual IR spectra. The strongest spectral lines are identified. Different models 
are indicated as in Fig. 2. 
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Fig. 5. — 4 — 2 0(3) line emissivity (cm~^ s"^ sr~^) as a function of the optical thickness 
(in mag) computed for model (i) (solid thick line) and model (vi) (solid thin line). The 
relative contributions to the population of the upper state of the transition, {v, J) = (4, 1), 
are: formation pumping (dotted line), IR ro- vibrational cascade within the ground electronic 
state (dashed line), UV fluorescent cascade from excited electronic levels (dot-dashed line), 
and thfirma.1 rnllisinns fdnt- dnt— dnt— dashftd linfiV Tinner na.npl" tra.nsliirftnt riniid mndfil ^ 
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Fig. 6. — 1 — S(7) line emissivities (cm~^ sr~^) as functions of the optical thickness 
(in mag) computed for models (ii), {iii), {iv — c), and {v — c) (solid lines), in the case 
of RCM. The relative contributions to the population of the upper state of the transition, 
{v,J) — (1,9), are: formation pumping (dotted line), IR ro-vibrational cascade within the 
ground electronic state (dashed line), UV fluorescent cascade from excited electronic levels 
(dot-dashed line), and thermal collisions (dot-dot-dot-dashed line). 
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Fig. 7. — Residual IR spectrum for a translucent cloud with gas density rin = 10^ cm ^ and 
Av = 5 mag. All other parameters are as for RCM. 
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Fig. 8. — Flux at Earth of the 4 — 2 0(3) transition as a function of H2 ionization rate for the 
RCM using model (i). Solid horizontal hues represent UKIRT sensitivity in the H (1.6 /xm) 
and K (2.2 //m) bands. 



